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We study how the Mott metal-insulator transition (MIT) 
is affected when we have to deal with electrons with different 
angular momentum quantum numbers. For that purpose we 
apply ab-initio quantum-chemical methods to lithium rings in 
order to investigate the analogue of a MIT. By changing the 
interatomic distance we analyse the character of the many- 
body wavefunction and discuss the importance of the s—p or- 
bital quasi-degeneracy within the metallic regime. The charge 
gap (ionization potential minus electron affinity) shows a min- 
imum and the static electric dipole polarizability has a pro- 
nounced maximum at a lattice constant where the character 
of the wavefunction changes from significant p to essentially s- 
type. In addition, we examine rings with bond alternation in 
order to answer the question under which conditions a Peierls 
distortion occurs. 



I. INTRODUCTION 

When the lattice constant of a metallic system is suf- 
ficiently enlarged it will become an insulator, before 
eventually we end up with a collection of well sepa- 
rated atoms. This metal to insulator transition was first 
pointed out by Mott([ 1] and references therein) and has 
stimulated a tremendous amount of work for a better 
understanding of it. The standard model for these in- 
vestigations has been the Hubbard Hamiltonian^ . It as- 
sumes one orbital per site and the electronic Coulomb 
interactions are reduced to an on-site Coulomb integral 
U. By applying different approximations Hubbard was 
able to quantify the conditions for the occurrence of a 
M-I phase transition. For that reason it is often called a 
Mott-Hubbard transition. More recently, the dynamical 
mean-field theory (DMFT) has been applied to the Hub- 
bard Hamiltonian with one electron per site and has fur- 
ther elucidated details of the transition (for an overview 
and further references see, e.g., [3]). 
The purpose of the present work is to go beyond the 
Hubbard model in order to understand better how other 
features might influence the transition in a realistic sys- 
tem. Lithium is a metal with a strong contribution of p 
electrons to the conduction band. With increasing lat- 
tice constant the ratio of p and s electrons is changing 
and in the atomic limit we are dealing with a pure ls^2s 
atomic state. One of the questions we would like to have 
answered is the effect of low-energy atomic excitations 
from 2s to 2p on the M-I transition. In order to come 



close to a realistic description of metallic systems we in- 
vestigate a ring of Li atoms applying quantum-chemical 
methods for a treatment of electronic correlations. 
In terms of a Hubbard model the Li-Li distance would 
correspond to the hopping term t. Instead of only an on- 
site Coulomb repulsion as in the Hubbard model we take 
here into account all possible Coulomb matrix elements. 
But in the limit of large seperations of atoms only the 
on-site Coulomb repulsion remains. 

The largest ring which we can treat at the multi- 
reference configuration-interaction (MRCI) level, with- 
out excessive computations, consists of 10 Li atoms. A 
medium-size basis set consisting of Gaussian type or- 
bitals (GTO's) is used in the calculations. Usually the 
onset of insulating behaviour is characterized by the ap- 
pearance of a gap. However, with a finite number of Li 
atoms in the ring we always have a finite energy dif- 
ference between consecutive energy levels. Therefore, 
when the atomic distance is changed we must look for 
a change in the pattern between the highest occupied 
molecular orbital (HOMO) and the lowest unoccupied 
orbital (LUMO). These quantities are only accessible 
at the mean-field level. At the correlated level, we se- 
lect three quantities which can give insight as regards 
the transition behaviour. MRCI methods can yield very 
good approximations to the exact ground-state wavefunc- 
tion. The character of this wave function will be the first 
criterion. A second quantity which we calculate is the 
one-particle energy gap defined as the energy difference 
which is obtained if an electron is added to, or removed 
from the system. As a third quantity, we determine the 
static dipole polarizability of the system by applying a 
small homogeneous electric field. These three quantities 
arc determined as functions of the nearest-neighbour Li- 
Li distance in the ring, and we discuss how the metal- 
insulator transition shows up in them. For these quasi- 
one-dimensional systems, the question arises whether we 
can gain energy by dimerization (Peierls distortion). We 
address this question by allowing for a bond alternation. 
In the next section we will discuss some technical details 
(Sect. II). In Sect. HI, we present the results of different 
mean-field methods, for Li rings with increasing number 
of Li atoms and with increasing Li — Li distance. The 
behaviour of the many-body wavefunction is discussed 
in Sect. IV. Gap energies are presented in Sect. V and 
polarizabilities in Sect. VI. The Peierls distortion is con- 
sidered in Sect. VII and conclusions follow in Sect. VIII. 
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Basis set 


IP (a.u.) 


EA (a.u.) 


a (a.u.) 




aug-cc-pVDZ 


0.196 


0.0151 


166.0 


2.73 


aug-cc-pVTZ 


0.196 


0.0176 


167.2 


2.69 


aug-cc-pVQZ 


0.196 


0.0209 


167.5 


2.68 


[4sl 

L "J 


0.196 


0.0019 


0.000 


3.03 


[3slp] 


0.196 


0.0150 


164.8 


2.85 


[4s Ip] 


0.196 


0.0212 


165.6 


2.85 


[4s2p] 


0.196 


0.0216 


165.6 


2.73 


[5s2p] 


0.196 


0.0218 


166.6 


2.67 


exp. 


0.198^' 


0.0227^' 


164.8'-' 


2.67^' 



TABLE I. ACPF values of the ionization potential (IP), 
the electron affinity (EA), and the dipole polarizibility (a) of 

the Li atom arc listed for different basis sets. Corresponding 
values for the bond length (ddimer) of the Li2 molecule are 
listed, too. 



II. TECHNICAL DETAILS 

The application of quantum-chemical ab-initio meth- 
ods to Li rings is at the focus of the present work. 
Starting from a mean-field Hartree-Fock (HF) descrip- 
tion, we re-optimize the valence wavefunction in multi- 
configuration sclf-consistcnt-ficld (MCSCF) calculations, 
keeping the Is core electrons frozen at the HF level. 
Thereby, the number of valence active-space orbitals is 
chosen equal to (or larger than) the number of Li atoms 
of the system, and all possible occupations of these or- 
bitals are accounted for. With this choice, we can prop- 
erly describe the dissociation limit where each atom has 
one valence electron. On top of the MCSCF calculations, 
we apply the multi-reference averaged coupled pair func- 
tional (MRACPF) method^'^. Here, aU configuration- 
state functions are included which can be generated from 
the MCSCF reference wavefunction by means of single 
and double excitations from the active orbitals. These 
calculations are performed with the program package 
MOLPRO^-". 

The GTO one-particle basis set used for the MC- 
SCF/MRACPF calculations of the present work is de- 
rived from the correlation-consistent polarized valence 
double-zeta (cc-pVDZ) basis set of Dunning and co- 
workers-'^^. From this set, we select the three s fimctions 
and the first contracted p-function. To describe the neg- 
atively charged system, we add an even-tempered diffuse 
s-function (exponential parameter: 0.0107) resulting in a 
[4s Ip] basis set. We checked the quality of this basis by 
calculating the electron affinity (EA), the ionization po- 
tential (IP) and the dipole polarizability of the Li atom 
(Table I). The p function is essential to describe the 
EA and the polarizability properly. Whereas the IP is 
weakly dependent on the size of the basis set used, be- 
cause correlation effects are absent (core- valence correla- 
tion is neglected), the EA substantially increases when a 
diffuse s-function is added. A second diffuse s-function 
or a diffuse p-function have nearly no influence. The po- 
larizability only slightly depends on the basis-set quality 



a = ao 


EA 


IP 


HF 


0.0208 


0.1554 


MCSCF 


-0.0027 


0.1672 


ACPF 


0.0335 


0.1812 


FCl 


0.0321 


0.1820 


a = L.'xio 


EA 


IP 


ACPF 


0539 


0.1530 


FCl 


0.0500 


0.1548 


a — 2.0ao 


EA 


IP 


ACPF 


0.0448 


0.1703 


FCl 


0.0409 


0.1707 



TABLE II. The EA and IP of Lig are listed for different 
levels of approximations and different Li-Li distances. 



beyond cc-pVDZ, if the basis set contains at least one p 
function. The errors of the selected [Aslp] basis set as 
compared with experimental data are 1% for the IP, 7% 
for the EA and 0.5% for the polarizability. Thus, the 
[4s Ip] basis although being relatively small can describe 
with sufficient accuracy all quantities we are interested 
in and is taken as default basis set in the subsequent cal- 
culations. 

The Li Li distance in 3-dimensional Li metal is 
ao=3.03A^, while the equilibrium distance of the Li2 
dimer (2.67A^'^) is smaller by about 10% . For the Liig 
ring, we calculate an equilibrium distance of 3.06A which 
is very similar to the bulk nearest-neighbour distance. 
Note, however that the [4slp] basis overestimates the Li2 
bond length by 0.18 A (the lack of a second p function 
is mainly responsible for this error, cf. Table I); if this 
bond- length error carries over to Liio, we would have an 
equilibrium distance of around 2.9 A for this system. 



III. MEAN-FIELD DESCRIPTION 

As a first step of our MIT study for Li rings, we per- 
formed mean-field (Hartrec-Fock and density-functional 
(DFT)) calculations for various ring lengths and internu- 
clear distances. With the Crystal program^^ we can do 
HF calculations for the one-dimensional infinite Li chain. 
As basis set we selected the optimized [3s2p] basis set 
for the three-dimensional metal^^. In Fig.l, the Hartree- 
Fock (HF) band structure is shown for the bulk equilib- 
rium internuclcar distance ao, and for a = 2.5ao which 
corresponds to the dissociation limit. In both cases the 
Fermi energy lies within the 2s band, but the band width 
is significantly reduced with increasing Li — Li distance. 
The predominant 2p band is very flat near the dissocia- 
tion limit and well separated from the 2s band, whereas 
for a distance ao there is no gap between the 2s and 2p 
band. 

For finite systems we applied the Turbomole ab-initio 
program package^^ to obtain mean-field results for very 
long rings up to n = 170. We explicitly imposed spatial- 
symmetry restrictions on our closed-shell wavefunctions 
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FIG. 1. Hartrec-Fock band structure for the infinite Li 
chain for two different Li — Li distances. The Fermi energy 
is set to zero. In addition, a band structure obtained from 
the levels of Liio is plotted (dashed curve); for comparison, 
the center of the 2s band is shifted to zero. 
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FIG. 2. Mean-field results (HF and LDA) for finite Li rings, 
in dependence of the number of atoms in the ring. The Li — Li 
distance is fixed at the metallic value (o = oo). 



FIG. 3. Mean-field results (HF in Dioh symmetry, 
DFT(LDA and GGA) in D-zh symmetry) for a Liio ring, in 
dependence of the Li — Li distance. The gap energy is taken 
as the difference of IP and EA; in the HF case, also the 
HOMO-LUMO gap value is given. 

in order to avoid symmetry breaking. (For rings with 
n >30, the HF approximation favours a symmetry bro- 
ken ground state.) If we generate a HF band structure 
from the orbital energies of a Liio ring and shift the Fermi 
energy to zero (dashed curve in Fig.l), the band width 
and band shape correspond well to that of the infinite 
chain, for both a = aa and a — 2. Sao- 
The energy gap of Li rings with different number of 
atoms, at a fixed Li — Li distance of a = ao, is shown 
in Fig. 2. In a mean-field approach, the gap energy can 
either be calculated as HOMO-LUMO gap in the neutral 
system or by adding and subtracting an electron to the 
system, i.e., as difference of the IP and EA. At the HF 
level both approaches (approximately) coincide accord- 
ing to Koopmans' theorem^''. Slight differences occur 
for smaller rings due to relaxation of the orbitals when 
the electron number is changed. This effect is missing 
in the HOMO-LUMO gap. Calculating the IP-EA gap 
with density functional theory in the local-density ap- 
proximation (LDA)^^, the results coincide well with the 
HF results. Note that in LDA the HOMO-LUMO gap is 
much smaller. This gap is directly related to an excita- 
tion energy other than the IP-EA gap. The decay of the 
IP-EA gap with increasing number of atoms in the ring 
is slow (n""'^^). (If wc would regard the finite ring as 
consisting of free electrons in a box, the decay should be 
proportional to (box- length)"^.) Despite the slow clos- 
ing of the gap, even for Liio it lies well within the band 
width of the infinite chain. 

For the Liio ring, we also calculated the behaviour of 
the IP-EA gap with increasing Li — Li distance, from 
ao to 2.5ao (Fig. 3). For the mean-field approaches, 
HF and DFT (both LDA^^ and GGAi^-^O), the gap is 
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monotonously closing. Thereby, the differences between 
the DFT variants and the differences between DFT and 
HF are not larger than that between the HOMO-LUMO 
and IP-EA values at the HF level. (Small differences may 
also arise due to the use of different formal symmetries 
(D2h for DFT, Dio/i for HF).) If we would model the ring 
by free electrons in a box again, the decay of the HOMO- 
LUMO gap should be proportional to (box-length)"^. 
This is in contrast to an approximate a^^ decrease ob- 
tained in our mean- field calculations, which shows that 
the electrons are far away from the free-electron limit. It 
should be kept in mind, however, that the correct dis- 
sociation limit cannot be described with these methods, 
anyway. Only multi-configuration methods can properly 
describe both the dissociation limit and the metallic be- 
haviour. Such calculations form the main part of our 
present work and will be discussed below. 



o ACPF, only s basis 
□ ACPF, metallic solution 
c ACPF, insulating solution 
y. ACPF, s start 



« □ 



« * * 



1.5 2 
Li-Li distance in multiples of a^ 



IV. THE CHARACTERISTIC FEATURES OF 
THE GROUND-STATE WAVEFUNCTION 

Taking the example of the Liio ring, we want to study 
how the character of the many-body ground-state wave- 
function changes when enlarging the Li — Li distance from 
a = ao to a = 2.5ao. For that purpose we performed 
a MCSCF calculation'''^ for the lowest singlet state, on 
top of a closed-shell HF calculation like that described 
in the previous section. (Using the MOLPRO package, 
we had to formally change to T>2h symmetry, however.) 
The selection of the active space is a decisive step: It 
is well known that for a correct description of the disso- 
ciation limit in Liio ten atomic 2s orbitals are needed, 
but for the metallic regime (a « ao) it is not clear how 
many orbitals are really important. Therefore we per- 
formed a pilot MCSCF calculation for a = ao with a 
high number of active orbitals and selected those natural 
orbitals whose occupancies were above a certain thresh- 
old. Retaining orbitals with occupancies larger than 
0.18 and rejecting all the others which have occupan- 
cies smaller than 0.03 resulted in 10 active orbitals, too. 
However, only the 9 energetically lowest ones fall into the 
same irreducible representations (in symmetry) as 
the 9 lowest-lying orbitals of the insulating regime, while 
the highest natural orbital of the metallic and insulating 
regime differ in symmetry. Therefore, we selected for our 
final MCSCF calculation a new set of 11 active orbitals 
which was generated as the union of the two sets of or- 
bitals described above. These orbitals were reoptimized 
at the MCSCF level. Finally, we performed a MRACPF 
calculation^'^ on top of the MCSCF, in order to include 
dynamical correlation effects as fully as possible. Note 
that the multi-reference treatment is important in both 
limits: in the metallic regime to describe static corre- 
lation with the quasi-degenerate 2p orbitals, and in the 
insulating regime to get the dissociation into degenerate 
2s orbitals right. 

Although we included in the active space all orbitals 



FIG. 4. MRACPF energies of the Liio ring versus the 
Li — Li distance. The curves labeled "insulating wavefunc- 
tion" and "metallic wavefunction" refer to different MCSCF 
zeroth-order wavefunctions used as reference for the subse- 
quent MRACPF, cf. text. In the inset an enlarged part of the 
transition region is shown. The filled symbols correspond to 
the MCSCF values. The curves labeled "s start" and "only 
s basis" refer to calculations leaving out p functions in the 
MCSCF zeroth-order wavefunctions and at both the MCSCF 
and MPACPF levels, respectively. 




1.4 1.6 
Li-Li distance in multiples of a„ 



FIG. 5. The p occupation of the Li atom in a Liio ring 
versus the Li-Li distance. The p occupations are calculated 
with the MuUiken population analysis. 
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which arc important in the hmits of small and large 
atomic distances, the total energy as function of the Li — 
Li distance is not a smooth curve (Fig. 4). This is so 
because the nature of the active orbitals changes along 
the curve. When starting the calculations from large a 
(wavefunction for the insulating state) and always us- 
ing the previous solution as starting point for the next 
smaller lattice constant, this yields a slightly different so- 
lution in the region of the MIT than when starting from 
the metallic regime and increasing a. The difference is 
larger, of course, for the MCSCF energies than at the 
MRACPF level; With a full CI (i.e., including all higher 
excitations) the differences would vanish. However, a sin- 
gle smooth curve can be obtained even at the MRACPF 
level if it is based on a MCSCF calculation with only s 
functions in the active space. In that case p excitations 
appear at the MRACPF level only. It is clear that such 
a procedure favours the insulating solution: it yields a 
ground state energy for a = ao which is by ~0.02 a.u. 
higher than the value obtained with variationally opti- 
mized MCSCF orbitals. Still, the so obtained ground- 
state potential curve is qualitatively correct also in the 
metallic regime. It is interesting to compare it with a 
MRACPF curve obtained by leaving out p functions in 
the one-particle basis set altogether. Here, we still ob- 
serve a (shallow) minimum near a = ao, but the energy 
is higher by more than 0.1 a.u. Thus, p contributions 
amount to about two thirds of the binding energy of the 
Liio ring. 

Turning back now to our most accurate potential curve: 
The change of the character of the MRACPF wavefunc- 
tion as well as that of the underlying MCSCF wave func- 
tion as a function of a is seen best by a Mulliken popula- 
tion analysis for the p orbitals (Fig. 5). For the "metal- 
lic" MCSCF solution the atomic p orbital occupancy is 
nearly 0.5 at a = ao and slowly drops to 0.3 at a w 1.5ao, 
while for the "insulating" solution it is less than 0.03 for 
a > 1.5ao. Even in the MRACPF case, the p population 
changes by a factor of ~3. Taking this rather dramatic 
change as a signature of the MIT we find for all investi- 
gated rings a region around 1.5ao , where the character 
of the wavefunction changes significantly. 



V. THE ONE-PARTICLE ENERGY GAP 

The energy gap of an insulator is determined by 

the energy differences between the groimd state of the 
neutral system and the ones with one electron added 
and subtracted. In the case of a Li ring, we calcu- 
late the MRACPF ground-state energies of the Li~ , Li„ 
and Li+ systems and determine EA=E(Li„)-E(Lijj^) and 
IP=E(Li+)-E(Li„) and from IP-EA the corresponding 
gap. For a metallic solid adding and removing an elec- 
tron cost an energy given by the chemical potential. In 
finite systems, there always remains the difference be- 
tween the one-particle energies of the additional and the 
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Li-Li distance in multiples of a^ 



FIG. 6. HF and MRACPF values for the EA and IP of a 
Liio ring, versus the Li-Li distance. For comparison, we add 
results from calculations without p functions in the basis set 
(" only s basis" ) . 



missing electron as well as the influence of relaxation and 

correlation effects. 

The EA and IP for Liio are plotted in Fig. 6. In the 
metallic regime, the HOMO-LUMO gap is expected to 

decrease with increasing Li-Li distance exactly like it 
happens for free electrons in a box with changing length. 
In fact, at the HF level we obtain a monotonous shrink- 
ing of the gap, both with our standard basis set ([4slp]) 
and with a pure s basis set ([4s]). The availability of p 
polarization functions stabilizes the charged systems as 
compared to the neutral one, leading to a slightly smaller 
gap for a = ao than with a pure s basis, but the dif- 
ference goes to zero for a ^ oo and the general trend 
of a closing gap with enlarging Li-Li nearest-neighbour 
distance is the same in both cases. However, the trend 
completely changes when correlation effects are included. 
The most important change concerns the a — > oo limit, 
where HF fails to describe correctly the dissociation. In- 
stead of a closing of the gap, the correct limit is the 
difference between the atomic IP and EA. This limit 
is obtained indeed, in our MRACPF calculations. As 
already noted in Sect. 2, the atomic EA is severely 
underestimated in the case of a pure s basis set, but 
this is efficiently remedied with our standard [4slp] basis 
set. We furthermore note that the EA converges much 
slower to the atomic limit than IP (EA(a = 2.35ao)=1.6 
EA(atom),IP(a = 2.35ao)=0.94 IP(atom)). We think 
that this is due to the larger size of the Li~ ion than of the 
Li"*" ion. Interestingly, the variation of IP and EA from 
a = ao to a — *■ 00 is not completely monotonous. Already 
for the case of a pure s basis, the gap is nearly constant 
(and even slightly decreases) in the range ao < fl < 1.3ao, 
and only then gradually opens up to the atomic limit. In- 
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FIG. 7. HF, MCSCF and MRACPF one particle energy 
gaps are plotted versus the number of Li atoms in the ring for 
two different lattice constants, cf. text. 

eluding p contributions leads to a signifieant increase of 
IP and EA for a = oq as well as to a distinctly non- 
monotonous behaviour of the gap as a function of a, 
with a minimum around a = 1.5ao. The effect is due 
to a stabilization of the charged systems, with respect to 
the neutral one, by means of polarization clouds around 
holcs/clcctrons. Thus, for a finite ring in the metallic 
regime the gap is closing when enlarging the lattice con- 
stant, at the MIT point it is minimal, then it starts to 
increase gradually approaching the atomic limit value. 
We investigated the gap behaviour not only for the Liio 
ring, but also for the Li2 dimer, for the Lig and Lig rings, 
and for the Lis cube. (Li4 arranged in a square has a 
negative EA; hence, the calculation of the one-particle 
energy gap does not make any sense. Li rings with odd 
number of Li atoms do not have a closed shell HF ground 
state and therefore are not considered here.) For Lii2, we 
performed calculations at a = ag and a — a^iT' where 
^MiT defined as the point where the character of the 
MCSCF wavefunction changes from significant p charac- 
ter to an essentially s one. For all investigated ring sys- 
tems, the minimal value of the gap occurs for nearly the 
same lattice constant in different rings ( between 1.42ao 
and 1.55oo). This distance coincides approximately with 
the one where the character of the wavefunction changes. 
In Fig. 7 the energy gap at a = ao and a = ajviiT is listed 
for different methods and for different ring lengths. When 
increasing the number of atoms in the ring, the minimal 
gap energy at the ACPF level decreases faster than the 
gap at a = ao. Unfortunately, the behaviour of the gap 
for much longer rings cannot be calculated by the meth- 
ods we use because it is computationally too expensive. 
The gap calculated at the MCSCF level seems to satu- 
rate for longer rings, i.e. only the dynamical correlations 



FIG. 8. The MRACPF static dipolo polarizability of the 
Li2 dimer, for the electric field parallel and perpendicular to 
the bond direction, is plotted versus the Li — Li distance. 

are closing the gap. As yet, we cannot make a reason- 
able extrapolation and find out for which ring length a 
zero gap is to be expected. But for all investigated rings 
we found similar behaviour with increasing Li-Li distance 
independent of the number of Li atoms in the ring. In 
spite of this we reemphasize that although we have a fi- 
nite system instead of a real metallic system, we still see 
a transition from a metallic-like regime (where the gap 
is closing with increasing Li-Li distance) to an insulat- 
ing regime where the gap is opening towards the atomic 
limit. 

In order to validate the methods used in the present pa- 
per, we performed fuU-CI (FCI) calculations^^ for Lig 
at three different lattice constants, i.e., we determined 
the lowest variational energies in each case which can 
be obtained with our standard [4slp] one-particle ba- 
sis set. The resulting EA and IP are listed in Table II. 
The ACPF values coincide well with the FCI results and 
therefore are expected to provide a reliable, quantitative 
description of the gap energy. 

Finally, let us note, that we observe a similar minimal gap 
for the Lig cube, where no sharp jump in the p occupa- 
tion pattern of the wave function is visible. In the cube, 
the minimal gap occurs at a somewhat larger distance 
(1.65ao) than it does in the ring and is also a little bit 
larger (0.09 a.u.), but the qualitative behaviour remains 
the same. 



VI. DIPOLE POLARIZABILITY 

As a third quantity which can be used to define the 
MIT, we calculate the static electric dipole polarizability 
of the system. We apply a static electric field of strength 
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FIG. 9. The static dipolc polarizability of the Liio ring, for 
the electric field in the ring plane, is plotted versus the Li — Li 
distance for two different methods. 



up to 0.035 a.u. and perform a quadratic fit for the en- 
ergy of the system, E{£) = E{0) — ^a£^, yielding the 
polarizability a. Higher order terms can be neglected 
for the small field strengths we chose. For a metallic sys- 
tem the polarizability should be infinite, therefore a steep 
increase of the polarizability can indicate an insulator- 
metal transition. 

As a test example, we selected the Li2 dimer and ap- 
plied the £ field both in the bond direction and in a 
direction perpendicular to that, while gradually enlarg- 
ing the Li — Li distance (Fig. 8). At the equilibrium bond 
length, we obtain an average polarizability of 118 a.u. 
(the experimental value is 111 a.u.^'*). If the electric is 
parallel to the bond direction, the polarizability exhibits 
a maximum at a Li — Li distance similar to that where the 
minimal gap occurs. For a field direction perpendicular 
to the bond, the polarizability increases monotonously 
to two times the atomic value as the HOMO-LUMO gap 
becomes smaller. 

For the Liio ring, we find the same behaviour as in the 
dimer when we apply the field perpendicular to the ring 
plane. When we apply the field in the ring plane we find 
a pronounced maximum. In Fig. 9, we plot the polariz- 
ability per atom for the Liio ring, both at the MCSCF 
and MRACPF levels. A crossing of the curves occurs 
within the transition region, since we have two differ- 
ent solutions at the transition. MCSCF and MRACPF 
polarizability values differ only slightly. At a = ag the 
calculated MRACPF value is 94.4 a.u.. Experimental 
data are available for the static dipole polarizability of a 
Liio cluster^'^; unfortunately the cluster geometry is not 
fully known, but it can be supposed to be two- or three- 
dimensional. The experimental value is 70.2 a.u. and 
it roughly agrees with our value. As already noted, the 
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FIG. 10. The total energy of the Liio ring is plotted versus 
the dimerization angle (defined as deviation from the ideal 36° 
angle in Liio, in degree) for different average lattice constants. 



polarizability increases with decreasing lattice constant 
(i.e., when coming from the insulating side). For an infi- 
nite system the polarizability should diverge at the MIT. 
We see a precursor of this divergence and can therefore 
use it to define a MIT even for a finite system. The rapid 
decay of the polarizability in both the metallic and the 
insulating regime is due to the opening of the gap. 



VII. LATTICE DISTORTION 

In this last section, we allow for a bond alternation 
along the ring, so that the Li atoms can form dimers. 
For average lattice constants larger than that where the 
M-I transition occurs (i.e., for lattice constants where the 
insulating solution is the ground state) the formation of 
dimers is energetically favoured. The resulting dimer in- 
ternuclear distance agrees well with the calculated dis- 
tance of the free dimer (2.88 A). The energy gain is due 
to bond formation between two Li atoms. In the metallic 
regime (a < aun) a Peierls distortion^^ occurs, stabiliz- 
ing the dimerized state (see Fig 10). However, decreasing 
the lattice constant further, we found a point (a = 1.3ao 
for Liio), below which the equidistant arrangement is the 
ground state, which would imply a metallic state for the 
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infinite chain. 

At the HF level, a bond alternation is found for all (av- 
erage) lattice constants; electronic correlations seem to 
suppress this bond alternation at the equilibrium lattice 
constant. This is in contrast to poly-acetylene, where 
correlations reduce the bond alternation without fully 
suppressing it^^. 



VIII. CONCLUSIONS 

We investigated the analogue of the metal-insulator 
transition for one-dimensional lithium, applying high- 
level quantum-chemical ab-initio methods. The aim of 
this investigation was to find out how the MIT is modi- 
fied from the one in the Hubbard model, when we have 
to deal with several orbitals per site. We have chosen Li 
as an example since here the different orbitals correspond 
to different angular momentum quantum numbers, i.e., s 
and p electrons. By treating a Liio ring we could apply 
accurate quantum-chemical ab initio methods. Following 
the original ideas of Mott, a metallic system will become 
an insulator when the lattice constant is sufficiently large. 
But in distinction to the Hubbard model it is here the 
interplay of the s and p orbitals which strongly affects 
the MIT. At the transition point (omit ~ 1.5ao) the 
character of the many-body wavefunction changes from 
significant p to essentially s character. Therefore it must 
be kept in mind that in a real solid the re-population 
of orbitals belonging to different angular momenta may 
have similarly strong influence on the MIT as changes in 
the ratio of the Hubbard U to the kinetic energy. We 
found that approximately at the same interatomic dis- 
tance where the p character of the wavefunction changes 
so strongly, the one-particle energy gap is minimal and 
the static electric polarizability has a maximum. Increas- 
ing the ring length, the minimal gap closes faster than 
the gap at the equilibrium distance; the reduced level 
spacing for both longer rings and larger lattice constants 
enhances dynamical correlation effects. Allowing for a 
bond alternation in the ring, the system dimcrizcs for 
all (average) lattice constants larger than ~ 1.3ao), for 
the insulating case due to Li2 molecule formation, for the 
metallic case due to Peierls distortion. For even smaller 
lattice constants, the correlations stabilize the equidis- 
tant arrangement of the Li atoms along the ring. 
Further investigations are planned for open boundary 
conditions (i.e., linear fragments of the infinite chain) 
and for other alkali metals (thereby, e.g. mixing different 
alkali atoms to mimic the ionic Hubbard model). 



APPENDIX A: LIST OF ABBREVIATIONS 

MIT metal-insulator transition 

HOMO highest occupied molecular orbital 

LUMO lowest unoccupied molecular orbital 

GTO Gaussian type orbital 

cc-pVDZ correlation consistent polarized valence double-zeta 

HF Hartree-Fock 

MCSCF multi-configuration self-consistent field approximation 

MRCI multi-reference configuration interaction 

MRACPF multi-reference averaged coupled pair functional 

IP ionization potential 

EA electron affinity 

DFT density functional theory 

LDA local density approximation 

GGA generalized gradient corrected approximation 

FCI full CI - configuration interaction with all excitations 
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